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Abstract 

Classes of multivariate and cone valued infinitely divisible Gamma distributions are 
introduced. Particular emphasis is put on the cone-valued case, due to the relevance 
of infinitely divisible distributions on the positive semi-definite matrices in applica- 
tions. The cone-valued class of generalised Gamma convolutions is studied. In partic- 
ular, a characterisation in terms of an Ito-Wiener integral with respect to an infinitely 
divisible random measure associated to the jumps of a Levy process is established. 

A new example of an infinitely divisible positive definite Gamma random matrix is 
introduced. It has properties which make it appealing for modelling under an infinite 
divisibility framework. An interesting relation of the moments of the Levy measure 
and the Wishart distribution is highlighted which we suppose to be important when 
considering the limiting distribution of the eigenvalues . 

Key words: infinite divisibility; random matrix; cone valued distribution; Levy process; matrix 
subordlnator. 
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1 Introduction 

The classical examples of multivariate and matrix Gamma distributions in the prob- 
ability and statistics literature are not necessarily infinitely divisible [14], l! 19ll . II40II . 
These examples are analogous to one-dimensional Gamma distributions and are ob- 
tained by a direct generalisation of the one-dimensional probability densities; see for 
example [15], 02311 . 11241. Working in the domain of Fourier transforms, some infinitely 
divisible matrix Gamma distributions have recently been considered in [5]|, II270 . Their 
Levy measures are direct generalisations of the one-dimensional Gamma distribution. 
The work of 027ll arose in the context of random matrix models relating classical and 
free infinitely divisible distributions. 

The study of infinitely divisible random elements in cones has been considered in 
O, II25II , II26I1 , II31I and references therein. They are important in the construction and 
modeling of cone increasing Levy processes. In the particular case of infinitely divisible 
positive-definite random matrices, their importance in applications has been recently 
highlighted in Q, ||8l, I128|| and (291. This is due to the fact that infinite divisibility 
allows modelling by matrix Levy and Ornstein-Uhlenbeck processes, which are in those 
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papers used to model the time dynamics of a dxd covariance matrix to obtain a so-called 
stochastic volatility model (for observed series of financial data). 

Generalized Gamma Convolutions (GGC) is a rich and interesting class of one-dimen- 
sional infinitely divisible distributions on the cone IR+ = [0, oo). It is the smallest class 
of infinitely divisible distributions on M+ that contains all Gamma distributions and that 
is closed under classical convolution and weak convergence. This class was introduced 
by O. Thorin in a series of papers and further studied by L. Bondesson in his book HIOII . 
The book of Steutel and Van Harn (39] contains also many results and examples about 
GGC. Several well known and important distributions on M+ are GGC. The recent survey 
paper by James, Roynette and Yor [ 16 ] contains a number of classical results and old 
and new examples of GGC. The multivariate case was considered in Barndorff-Nielsen, 
Maejima and Sato [SJ. 

There are three main purposes in this paper. We formulate and study multivariate 
and cone valued Gamma distributions which are infinitely divisible. Second, we con- 
sider and characterise the corresponding class GGC{K) of Generalised Gamma Con- 
volutions on a finite dimensional cone K. Finally, we introduce a new example of a 
positive definite random matrix with infinitely divisible Gamma distribution and with 
explicit Levy measure. 

The main results and organisation of the paper are as follows. Section 2 briefly 
presents preliminaries on notation and results about one-dimensional GGC on M+ as 
well as some matrix notation. Section[3]introduces a class of infinitely divisible d-variate 
Gamma distributions Tj_{a, (3), whose Levy measures are analogous to the Levy measure 
of the one-dimensional Gamma distribution. The parameters a and /3 are measures and 
functions on S (the unit sphere with respect to a prescribed norm), respectively. It is 
shown that the distribution does not depend on the particular norm under consider- 
ation. The characteristic function is derived and it is shown that the Fourier-Laplace 
transform on exists if /3 is bounded away from zero a— almost everywhere. Further- 
more, the finiteness of moments of all orders is studied and some interesting examples 
exhibiting essential differences to univariate Gamma distributions are given. 

Section|4]considers cone valued Gamma distributions and their corresponding class 
GGC{K) of Generalised Gamma Convolutions on a cone K, defined as the smallest class 
of distributions on K which is closed under convolution and weak convergence and con- 
tains all the so-called elementary Gamma variables in K (and also all Gamma random 
variables in K in our new definition). This class is characterised as the stochastic inte- 
gral of a non-random function with respect to the Poisson random measure of the jumps 
of a Gamma Levy process on the cone. This is a new representation in the multivariate 
case extending the Wiener-Gamma integral characterization of one-dimensional GGC 
on E+ — [0, oo), as considered, for example, in II16I1 . 

Section [5] considers the special cone valued case of infinitely divisible positive- 
semidefinite dxd matrix Gamma distributions. New examples are introduced via an 
explicit form of their Levy measure. They include as particular cases the examples con- 
sidered in O, 02711 . A detailed study is done of the new two parameter positive definite 
matrix distribution AT{ri, E), where 77 > (d — l)/2 and S is a d x d positive definite ma- 
trix. This special infinitely divisible Gamma matrix distribution has several modeling 
features similar to the classical (but non-infinitely divisible) matrix Gamma distribution 
defined through a density in particular the Wishart distribution. Namely, moments of all 
orders exist, the matrix mean is proportional to S and the matrix of covariances equals 
the second moment of the Wishart distribution. When S is the dxd identity matrix 
Id, the distribution is invariant under orthogonal conjugations and the trace of a ran- 
dom matrix M with distribution AT{ri, Id) has a one-dimensional Gamma distribution. 
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A relation of the moments of the Marchenko-Pastur distribution with the asymptotic 
moments of the Levy measure is exhibited. Hence, this matrix Gamma distribution has 
a special role when dealing with a random covariance matrix and its time dynamics, 
e.g. by specifying it as a matrix Levy or Ornstein-Uhlenbeck process. As an application, 
the matrix Normal-Gamma distribution is introduced, which is a matrix extension of the 
one-dimensional variance Gamma distribution of 11220 which is popular in finance. 

2 Preliminaries 

For the general background in infinitely divisible distributions and Levy processes 
we refer to the standard references, e.g. II36II . 

2.1 One-dimensional GGC 

A positive random variable Y with law /x = C{Y) belongs to the class of Generalised 
Gamma Convolutions (GGC) on K+ = [0,cx)), denoted by T(M_|_), if and only if there 
exists a positive Radon measure on (0, oo) and a > such that its Laplace transform 
is given by: 

Lf,{z) =Ee-^^ =exp(^~az~ In (^1 + ^) u^(ds) j (2.1) 

with 

|loga;|u^(da;) <C50, / ^^^^^ < C50. (2.2) 

For convenience we shall work without the translation term, i.e. with a — Q. The mea- 
sure is called the Thorin measure of /i. Its Levy measure is concentrated on (0, oo) 
and is such that: 

v^j^l&x) = x^'^l^{x)dx, (2.3) 
where is a completely monotone function in a; > given by 

/•oo 

Z,,(da;) = / e-^'v^{As). (2.4) 

The class T(M+) can be characterized by Wiener-Gamma representations. Specif- 
ically, a positive random variable Y belongs to r(M+) if and only if there is a Borel 
function h : M+ — > M+ with 

poo 

\n{l + h{t))dt<oo, (2.5) 
such that Y = has the Wiener-Gamma integral representation 

/■oo 

Y'' ^ / h{u)d-fu, (2.6) 
Jo 

where {'jt]t > 0) is the standard Gamma process with Levy measure v^dx) = e~^^. The 
relation between the Thorin function h and the Thorin measure v^j, is as follows: is 
the image of the Lebesgue measure on (0,oo) under the application : s — > l/h{s). That 
is, 

/•oo poo 

T^ds^ e-^^'u^(dz), 2:>0. (2.7) 
Jo 

On the other hand, if F^^^ (x) — w^(d?/) for a; > and F^J-{s) is the the right continuous 
generalised inverse of F^j^is), that is F~^{s) = inf{t > 0;Fv^{t) > s} for s > 0, then, 
h{s) = 1/ F-\s) for s>0. 
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Many well known distributions belong to r(M+). The positive a-stable distributions, 
< a < 1, are GGC with h{s) = {ser{a + for a > 0. In particular, for the 

1/2— stable distribution, h{s) — 4 (s^tt) ^ . Beta distribution of the second kind, lognor- 
mal and Pareto are also GGC, see flF]. 

For more details on univariate GGCs we refer to IIIOI I16II 

2.2 Notation 

Mrf(M) is the linear space oi d x d matrices with real entries and its subspace 
of symmetric matrices. By §J and §^ we denote the open (in S^) and closed cones of 
positive and nonnegative definite matrices in M(i(M). S^d ||.|| is the unit sphere on K'^ 
with respect to the norm || • ||. 

The Fourier transform /i of a measure /i on M = M'' or M = Mrf(M) is given by 



where we use {A,B) = tT{A^ B) as the scalar product in the matrix case, where 
denotes the transposed on Md(M). By we denote the d x d identity matrix and by |A| 
the determinant of a square matrix A. For a matrix A in the linear group QCd{R) we 
write A-^ = {A'^y^ ■ 

We say that the distribution of a symmetric random dxd matrix AI is invariant under 
orthogonal conjugations if the distribution of OMO^ equals the distribution of M for 
any non-random matrix O in the orthogonal group 0{d). Note that M — > OMO^ with 
O <E 0{d) are all linear orthogonal maps on (or M^^) preserving §j". 

3 Multivariate Gamma Distributions 

3.1 Definition 

Definition 3.1. Let /i jbe an infinitely divisible probability distribution on M''. If there 
exists a finite measure a on the unit sphere Sgd ||.|| with respect to the norm \ \-\\ equipped 
with the Borel a -algebra and a Borel-measurable function {5 : Sgd ||.|| — > such that 

/t(z) = exp I / / f e*™ " - Ij dra(dw) (3.1) 

for all z e M.^, then fi is called a rf-dimensional Gamma distribution with parameters a 
and 13, abbreviated Td{a, /3) -distribution. 

If /? is constant, we call ^ a || • || -homogeneous Td{a, /3)-distribution. 

Observe that the notation Td{a,P) implicitly also specifies which norm we use, be- 
cause a is a measure on the unit sphere with respect to the norm employed and (3 is 
a function on it. The parameters a and (3 play a comparable role as shape and scale 
parameters as in the usual positive univariate case. 

Remaric 3.2. (i) Obviously the Levy measure of is given by 

v^.{E) - / / lE{rv) Ara{dv) (3.2) 

for all E € B{R'^). This expression is equivalent to 

^-f3{x/\\xmx\\ 

i^f^idx) = —r 5(da;), a; G M'' (3.3) 

11-^11 
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where a is a measure on W'- given by 

a{E)= I I lE{rv)Ara{<lv), E e B{R'^). (3.4) 

(ii) Likewise we define M(i(M) and S^-valued Gamma distributions with parameters 
a and (3 (abbreviated rMa{c^,P) and Ts^{a,l3), respectively) by replacing M'' with Md(M) 
and Sd, respectively and the Euclidean scalar product with {Z,X) =tr(X^Z). All up- 
coming results immediately generalise to this matrix-variate setting. We provide further 
details in Section[5] 

If d = 1 and a({ — 1}) — 0, then we have the usual one-dimensional r{a{{l}), /3{1))- 
distribution. In general it is elementary to see that for d = 1 a random variable X ^ 
ri(a,/3) if and onlyifX = X1-X2 with Xi - r(a({l}), /3(1)) and - r(a({-l}, /3(-l)) 
being two independent usual Gamma random variables, i.e. X has a bilateral Gamma 
distribution as analysed in II17I I18II and introduced in QUI 12211 under the name vari- 
ance Gamma distribution. If a({l}) = Q!({— 1}) and /3(1) = /?(— 1), it indeed can be 
represented as the variance mixture of a normal random variable with an independent 
positive Gamma one (a comprehensive summary of this case can be found in L39.I where 
it is called sym-Gamma distribution). 

Now we address the question of which a, f3 we can take to obtain a Gamma distribu- 
tion. 

Proposition 3.3. Let a be a finite measure on SRdj|.j| and /3 : Sgd ||.|| M+ a measur- 
able function. Then ( |3.2| l defines a Levy measure and thus there exists a rd{a,l3) 
probability distribution ^ if and only if 

In M + U(dv) < 00. (3.5) 

~ ,1111 ^ P\'")/ 

Moreover, J^ai\\x\\ A l)i^^(dx) < 00 holds true. 

The condition ( |3.5[ l is trivially satisfied, if /3 is bounded away from zero a-almost 
everywhere. 



Proof. 



1 /"I — 
x\\i^Jdx) = / / e-''(")''dm(di;) = / — a(dw) 



a;||<l 



< a(S]{dj|.||) < 00 

using the elementary inequality 1 — e^^ < a;, for each x e R+. Denoting by Ei the 
exponential integral function given by Ei{z) = for z e M+, we get 



v^,{dx) = / / dra(di;) = / 



Ei{P{v))a{dv) (3.6) 
Ei{z)T{dz), (3.7) 

where we made the substitution z = j3{v) and t{E) a{j3^^{E)) for all Borel sets E in 
M+. Since r is a finite measure and < Ei{z) < e"^ ln(l + l/z) Vz e M+ (see IT] p. 229]), 

00 

Ei{z)T{dz) < 00. 

1/2 
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The series representation Ei{z) = —7 — ln(z) — J2iLi ^ with 7 being the Euler- 
Mascheroni constant (LL P- 229]) implies that lim^^o ln(2;)) = 1. Hence, 

.1/2 .1/2 .1/2 

/ Ex{z)T{dz) < CO \ I ln(z)|T(dz) < 00 / ln(l + l/2;)r(dz) < 00 
JO ^0 Jo 

using ln(l + 1/z) = ln(l + — ln(z) and the finiteness of r in the second equivalence. 
Appealing to the finiteness of t once more, the above conditions are equivalent to 

ln(l + l/z)T(dz) = / ln(l + l//3(w))a(dw) < CX3. 

K"*, 11-11 

□ 



"'S, 



The next proposition shows that the definition of a Gamma distribution does not 
depend on the norm, only the parametrisation changes when using different norms. 

Proposition 3.4. Let || • ||a Jbe a norm on and ^be a Td{a, (3) distribution with a being 
a finite measure on Sgd ||.||^ and (3 : SRd.||.||^ M+ measurable. If\\ ■ \\b is another norm 
on M'^, then fi is a Td{ab,Pb) distribution with ab being a finite measure on Sjjd ||.||j^ and 
Pb ■ Sgd ii ii^ M+ measuraWe. Moreover, it holds that 

ab{E)^ f 1e ( jr^] a{dv) V E e B{S^an.\L) (3.8) 



Mvb) = 13 (tT^) \\vb\\a yVb e SRd,,,.,,^. (3.9) 

\\\Vb\\a J 

The above formulae show that the mass in the different directions, which is given 
by a, does not change, and (3 only needs to be adapted for the scale changes implied by 
the change of the norm. 

Proof. Substituting first Vb — v/\\v\\b and then s = r/\\vb\\a gives: 

-I3(v)r 

/ e'™ ' ^ - 1 j - 



exp 



(e'™"^^ - 1) ^ — dra(di;) 

(^e'Wu"''' - 1) - — ^'^"'""^ drabidvb) 




1 dsabidvb) 



□ 



3.2 Properties 

In this section we study several fundamental properties of our Gamma distributions. 
Proposition 3.5. Any Td{a, (3) -distribution is self-decomposable. 

Proof. This follows immediately from the definition and [36, Th. 15.10]. □ 

Later on we will considerably improve this result by showing that we are in a very 
special subset of the self-decomposable distributions. This result has important impli- 
cations for applications where one likes to work with distributions having densities, i.e. 
distributions which are absolutely continuous (with respect to the Lebesgue measure). 
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Proposition 3.6. Assume tiiatsuppa is of full dimension, i.e. that it contains d linearly- 
independent vectors in K"^. Then the Td{a, 13) -distribution is absolutely continuous. 

Proof. It is immediate that the support of Tdia, /3) is the closed convex cone generated 
by suppa. Hence, the support of Td{a,/3) is of full dimension and so the distribution is 



It follows along the same lines that in the degenerate case the rd(a, /3)-distribution 
is absolutely continuous with respect to the Lebesgue measure on the subspace gen- 
erated by suppa. If suppa consists of exactly d linearly independent vectors, rd(a,/3) 
equals the distribution of a linear transformation of a vector of d independent univariate 
Gamma random variables with appropriate parameters and thus the density can be cal- 
culated easily using the density transformation theorem with an invertible linear map. 
If supp a is a finite set of full dimension, one can calculate the density from the density of 
independent univariate Gamma random variables by using the density transformation 
theorem with an invertible linear map and integrating out the non-relevant dimensions. 
In general the density can be determined via solving a partial integro-differential equa- 
tion (see [37]). Moreover, criteria for qualitative properties of the density like continuity 
and continuous differentiability can be deduced from the results of II33I I34II , but look- 
ing at the simple case of a vector of independent univariate Gamma distributions one 
immediately sees that the sufficient conditions given there are far from being sharp. 
Therefore we refrain from giving more details. 

Next we show that our d-dimensional Gamma distribution has the same closedness 
properties regarding scaling and convolution as the usual univariate one. 

Proposition 3.7. (i) Let X - Td{a, P) and c> 0. Then cX - Td{a, (3/c). 

(ii) Let Xi ^ Td{ai,l3) and X2 ^ Td{a2,P) be two independent d-dimensional Gamma 
variables. Then Xi + X2 ^ ^diai + a2,/3). 

Proof. Follows immediately from considering the characteristic functions. □ 

Likewise it is immediate to see the following distributional properties of the induced 
Levy process. 

Proposition 3.8. Let L be a Td{a,(3) Levy process, i.e. Li ^ Td{a,(3). Then Lt ^ 
Td{ta,P) for all te M+. 

Of high importance for applications is that the class of Td distributions is invariant 
under invertible linear transformations. 

Proposition 3.9. Let X ^ Td{a, (3) (with respect to the norm \\-\\) and A be an invertible 
d X d matrix. Then AX ^ rd{aA,l3A) with respect to the norm \\ ■ \\a = ||^^^ • || and 



non-degenerate. Thus [35] concludes. 



□ 





(3.11) 



(3.10) 



Proof. We have for all z e M' 



d 




,i<z,AX> 



) 





dra(dv) 



f,-l3A{n)r 



draA{du) 



r 
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where we substituted u — Av. □ 

It is easy to see that the above proposition can be extended to m x d matrices of 
full rank with m > d. Obviously, such a result cannot hold in general for a linear 
transformation A with ker(A) ^ {0}, since combinations of one dimensional Gamma 
distributions are in general not univariate Gamma distributions. 

Next we present an alternative representation of the characteristic function. 

Proposition 3.10. Let n be Td{a,l3) distributed. Then the characteristic function is 
given by 




/i(z)=cxp(/ In ( ^^^P ; J a{dv) ] forallzeR'^ (3.12) 



where In is the main branch of the complex logarithm. 
Proof. Follows from the definition and the well known fact 



e 



/3(w) — iv^ z 



□ 



Note that if a has countable support {vj}jGN/ then 



= n ( 



P {vj ) — iv- z 



We now show that the Fourier-Laplace transform of a Gamma distribution exists if 
and to a certain extent only if (5 is bounded away from zero a almost everywhere. 

Theorem 3.11. (i) The Fourier-Laplace transform jl of a Td{a,l3) distribution ^ exists 
for all z in a neighborhood U of zero, if (3{v) > k for v € S^d ||.|| a-a.e. with k > 0. 
fi is analytic there and given by formula l |3.12| l. 

(ii) If there exists a sequence (wn)neN in Sr^ ||.|| with lim„^oo P(vn) = and a{{vn}) > 
for all n £N, then the Fourier-Laplace transform p, exists in no neighborhood U (^C^ 
of zero. 



Proof. Using Proposition 3.4 we can assume w.l.o.g. that the Euclidean norm || • II2 is 
used for the definition of the Td{a, (3) distribution. 

(i) We will no w show (i) for U = BJO) C C^, where B«(0) := {x e C : |jx||2 < k}. 
From Proposition 3.10 it is clear that fi{z) exists for all z e -Bk(O) C C'', if and only if 

exists for all z e -8^(0). Consider now an arbitrary S e (0,1) and z e ^^^(O). Then 
the Cauchy-Schwarz inequality implies < ||2:||2 < 6k and hence < S. 

Therefore In ^1 — -^5^^ exists and is bounded on Bs^l^O) a-a.e. This implies that 

exists on Bst^{0). Since S S (0, 1) was arbitrary, this concludes the proof of (i), since the 
analyticity follows immediately from the appendix of [|12ll . 
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(ii) W.l.o.g. assume f3{vn) < 1/n. For n e N set z„ — — i/3(u„)u„. Then ||2;„||2 = I3{vn) < 
l/n and 1 - (iw^z„)//3(w„) — 0. Hence, 



In 1 - ——^ a(dv) 



and thereby 



do not exist. This implies that fi is not defined on i3i/„(0). Since n e N was arbitrary, 
this shows (ii). □ 

Remark 3.12. To extend this result to the matrix case, one simply has to use the Frobe- 
nius or trace norms and the scalar product Z,X i-^ iv{X^ Z) instead of the Euclidean 
norm and scalar product. We consider this in Section [5] 

Proposition 3.13. A Td{a,j3) distribution fj, has a finite moment of order k > 0, i.e. 

/jjrf ||a;||*''/i(da;) < oo, if and only if 

P{v)-''a{dv) < oo. (3.13) 

"''■ii-ii 

Moreover, if m is the mean vector and S — (cri^)i,j=i,...,ci is the covariance matrix of 

rd(«,/?) 

Pivy^vaidv). (3.14) 



and 

S = / ^{vy^vv'^a{dv) (3.15) 



Proof. If f3 is bounded away from zero, l |3.13| l holds trivially and Theorem |3.11| implies 
that ^ has finite moments of all orders fc > 0. So w.l.o.g. assume that (3 is not bounded 
away from zero in the following. By [ 36. p. 162] fi has a finite moment of order k, if and 
only if 



r*^ dra(dv) < oo. 



Substituting s = rf3{v) this is equivalent to 



/ s'''^e"'dsa{dv) < oo. (3.16) 
.11-11 •^^(-) 

Assuming without loss of generality that (3{v) < 1 for all v e S^d ||.||, we have that 



< C(fc) := 



nOO nOO 

/ s''-^e-'ds< / s'^^ie-Ms < r(fc). 

Jl J 



Hence, ( |3.16[ l is equivalent to l |3.13| l. Finally, l |3.14[ l and (jSTTSJ follow from Example 
25.12 in II36II and observing that that the infinitely divisible distribution Td{a,j3) with 
Fourier transform (|3.l|l has Levy triplet (C, 0, v^, where Q = xi'^{dx). □ 



Corollary 3.14. A || • \\-homogeneous Td{a,j3) distribution has an analytic Fourier-La- 
place transform in i?^(0) and finite moments of all orders. 
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Hence, any homogeneous Gamma distribution behaves like one would expect it from 
the univariate case. However, the behaviour in the non-homogeneous case may be 
drastically different, as the following examples illustrate. 

Example 3.15. Consider d~2. Let a be concentrated on {w„}„eN with 

Vn = (sin(n^"'"), cos(n~"'")) 

and set a{{v.n}) = e^" and /3(i'„) — 1/n for all n e N. Then by Theorem |3. 11 1 (ii) the 
Fourier-Laplace transform exists in no neighbourhood of zero. 

r« 0(v)^^a{dv) = CM n^e^^^ is finite for all fc > using the quotient criterion, 
because 

(n-l- l)'=e-("+i' , 
lim = e < 1. 

Thus, we have moments of all orders, but the Fourier-Laplace transform exists in no 
complex neighbourhood of zero. 

Example 3.16. Consider the set-up of Example 3.15 but set now a({w„}) — l/n^^™ for 



some real m > 0. /g ^ I3{v) ^a{dv) = X^neN ^i^^ finite if and only if A; < m. 

It is easy to see that condition ( |3.5| l is satisfied if condition ( |3.13| l holds for some k > 
0. Hence, the r2(a,/3) distribution exists indeed, but only moments of orders smaller 
than m are finite. 



Example 3.17. Consider again the set-up of Example 3.15 Set now a{{vn}) = (ln(l + 
nf{n+l))-'^. 

Then /^^^^ In (l + ^) a{Av) ^ Y.nen in(i+n)^(i+n) < °« (see (321 Theorem 3.29] 
and thus the T2{a, /3) distribution is well-defined. 

Yet, /s^, ^1 „^ /3(w)-'=a(dz;) = E„eN in(i+nt^(i+«) = °° fo^" ^11 real fc > and so the 
r2(a, /?) distribution has no finite moments of positive orders at all. 



4 Gamma and Generalised Gamma Convolutions on Cones 

4. 1 Cone-valued infinitely divisible random elements 

We first review several facts about infinitely divisible elements with values in a cone 
of a finite dimensional Euclidean space B with norm || || and inner product (•, ) . A 
nonempty convex set if of i? is said to be a cone if A > and x £ K imply Xx e K. A 
cone is proper if a; = whenever x and —x are in K. The dual cone K' of K is defined 
as K' = {y <E B' : {y, s) > for every s e K} . A proper cone K induces a partial order 
on B by defining xi <k x2 whenever X2 — xi <E K for xi e B and X2 G B. Examples of 
proper cones are M+, E^J. = [0,oo)'*, §J and 

A random element X in K is infinitely divisible (ID) if and only if for each integer 
p> 1 there exist p independent identically distributed random elements Xi, ...,Xp in K 

such that X ''= Xi + ... + Xp. A probability measure /Li on iiT is ID if it is the distribution 
of an ID element in K. It is known (see |381) that such a distribution /i is concentrated 
on a cone K if and only if its Laplace transform L^(9) = exp(— (9, x))iJ,{dx) is given 
by the regular Levy-Khintchine representation 

L^(e)=exp|-(e,*o>-y^(l-e-<®'">)j^^(da:)| for all 9 G if', (4.1) 

where G if and the Levy measure is such that Vp{K'^) = and 

/ (||x|| Al)i/^(da;) <oo. (4.2) 



Infinitely Divisible Multivariate Gamma Distributions 



11 



If X = {X{t);t > 0} is the _ftr-increasing Levy process (/-sT -valued subordinator) asso- 
ciated to n, its Levy-Ito decomposition is of the form 



X{t) = + / / : 
JQ JK 



xN{At,Ax) 



s<t 



where AX(s) e K for all s > a.s. and N{dt,dx) is a Poisson random measure on 

E{Af (dt, Ax)) = v^{Ax)At. (4.4) 



X K with 



4.2 Cone-valued Gamma distributions 

Let Sjj.ji be the unit sphere of B with respect to the norm || || and let X be a proper 
cone of B. We write S|^|| = S||.|| n K and denote by B{Sf^^ ) the Borel sets of Sl":'!!, 
Similar to Definition [STT] we have Gamma distributions in the cone K. 



Definition 4.1. Let ^ be an infinitely divisible distribution on the cone K. If there exist 
a finite measure a on Sj|^|| , and a measurable function (3 : Sjj^jj such that 

- I (l - e-'-^^'^)) dra{dU) \ (4.5) 

II " ^ ) 

for all Q £ K' , then fi is called a K -Gamma distribution with parameters a and (3, and 
we write /i ^ F/f (a, /3). The Levy measure of ^ is 

i^^{E)^ / lE{rU) dra{dU), E&B{K), (4.6) 

7s-„ Jo 

and satisfies 

min(l, |lx||)i^^(da;) < CX3. (4.7) 



IK 

The expression (|4.6|l is equivalent to 



-^)(X/||X||)||A'|| 

i^^idX) = — lKiX)aidX), (4.8) 



where a is a measure on K given by 



a{E)^ [ [ lE{rU)dra{dU), E e B{K). (4.9) 
•^^^11 ^0 

All properties of the multivariate Gamma distribution in Section [3] are also true for 
the cone-valued Gamma distribution. As in Proposition 3.3 we can in particular show 
that there exists a Txia, P) probability measure /_* if and only if 



In (^1 + ^ j a(dC/) < oo, (4.10) 

I - 1 

in which case we have l |4.7[ l. Also, as for Proposition 3.10 the Laplace transform of a 
Fif (a, j3) probability measure ^ is also given by 

L^(e)=exp|-^^ ln(^l + ^|^^a(dC/)|, 9 G X'. (4.11) 
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If II 'll J is another norm on K and if ^ has distribution (a, (3), then /i has distribution 
Tif /?;,) where a;,, are given as in ( |3.8| l and l |3.9| l respectively. Also, ri<:(a, /3), has a 
finite moment of order fc > 0, if and only if 

l3{U)-^a{dU) < oo. (4.12) 



In the homogeneous case, i.e. I3{U) — (3o > for any U e S|''|| , we have E || A/|| < oo for 
any fc > 0. 



If M is a random element in K with distribution TKia, /3) and (4.12 1 is satisfied with 

k = l, 

E{M) = [ l3{U)-^Ua{AU). (4.13) 

4.3 Ito-Wiener-Gamma integrals 

In this section we formulate an Ito-Wiener-Gamma integral for -valued Gamma pro- 



cess, similar to the Ito-Wiener-Gamma integral 1 2.6 1 with respect to the one-dimensional 
Gamma process. 

Let 7 — 7(0;, /3) = (7t;i>0) be a X-valued Gamma process. That is, 7 is the 
JsT-increasing Levy process such that ^aji = rif(a, /3) is the distribution of 71. Let 
N~f{ds,dx) be the random measure on M+ x K associated to the ii' -valued jumps of 
7 and v^^ j, be the Levy measure of 71. Hence E{iV(di, dx)} — i/^^^ (da;)dt, where 

r -fi{U)r 
^i^^fAE)^ / lE{rU) dra{dU), EeB{K). 

Let h : M+ x S^^,, Mi be a measurable function such that 



00 







In ( 1 + ^'^^ ) a(d;7)ds < cx), (4.14) 



in which case we say that h belongs to L{T k{(^, P))- The last condition is the cone 
analogon of the one-dimensional condition ( |2.5[ l. 

We prove in the next proposition that the following Ito-Wiener-Gamma integral type 
is well defined 

^""i ^/i(s,^)a;7V(ds,dx), (4.15) 

in the the framework of integration with respect to infinitely divisible independently 
scattered random measures (i.d.i.s.rm.) in Rajput and Rosinski L30.I (see also [9] for the 
special case of random matrices). 

Proposition 4.2. The integral l |4.15[ l is weJ] defined if and only if the function h : M+ x 
S^ll M+ belongs to L(rif (a, /?)). 

Moreover, h e L(ri<-(a, /?)) if and only if 

min(l, \\h{s, xj ||x||)a;||)t'^t^ ^ (dx)ds < cx), (4.16) 

JK 

and (4. 1 6p is equivalent to the following two conditions 

»l/2 

|hi(z)|G'[/(dz)a(d[/) < CX) (4.17) 

si^„ Jo 



Infinitely Divisible Multivariate Gamma Distributions 



13 



and ^ 

/ / -Gu{dz)a{dU) < oo (4.18) 

7sf „ Jl/2 Z 

where Gu{dz) is the measure on which is the image of the Lebesgue measure on M+ 
under the change of variable s — > j3{U)/h{s, U). 

Proof. By [|9l l30ll the existence of the integral is equivalent to l |4.16[ l. 

Since /3(t/) < a.e. U, Fubini's theorem and elementary computations give 



oo 



/ = / / min(l,||/i(s,a;/||a:|l)x|l):/^„,(d:E)ds (4.19) 

Jo JK 



oo /-oo „—rl3{U) 



min(l, ||/i(s, U)rU\\)- dm(dC/)ds 



"'0 

oo pl/h{s,U) 



"'0 



oo ^oo ^-rl3{U) 



s^ii Jo Ji/h{s,u) r 

m) 



draidU) 



1 - ) dsa{dU) 



Sm „ "'0 



El ( ] dsa{dU) (4.20) 



where Ei is the exponential integral function as in the proof of Proposition |3.3[ 
Using the change of variable z — j3{U)/h{s, U) we have 



- — —Gu{dz)a{dU) 



+ / / Ei[z)Gu{dz)dsa{dU) = Ii+ h (4.21) 

We shall in the following show that / < oo if and only if 

13=/" / ln(l + -jG[/(dz)a(df/)<oo (4.22) 

if and only if l |4.17[ l and ( |4.18[ l are satisfied. This concludes, as obviously ( |4.22[ l and 

h G L{TK{a, P)) are equivalent. 
First, 

r r°° 1- e""" 

/ / Gu{dz)a{dU) < oo (4.23) 

if and only if ( |4.18[ l, if and only if 



h= f I ln(l + -)G'[/(dz)a(dC/)<oo 

^Sf „ Jl/2 V Z J 



since (l/z)/ln(l + z ^) — > 1 as z — > oo. 
On the other hand, 

.1/2 



s f.. Jo 



El (z) Guidz)aidU) < oo 
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if and only if (4.17 1 holds (since Ei{z)/{— ln(z)) 1 as z 0) if and only if 

.1/2 / ^ 



s f„ Jo 



In 1 + - Gu{dz)a{dU) < oo 



because ln(l + z ^)/{— ln(z)) — > 1 as z — > 0. 
/s < oo and ( |4.17| l both imply 



sf „ JO 



1/2 



Gu{dz)a{dU) < oo. 



Thus 



s^„ Jo 



1/2 



1 - e- 



-Guidz)aidU) < 



sr. "'0 



1/2 



Guidz)a{dU) < oo, 



provided /s < oo or ( |4.17[ l hold. 



Since < Ei{z) < e"^ ln(l + 1/z) Vz e M+ (see IH p. 229]), using ( |4.23[ l implies 

El (z) G,7(dz)dsa(dC/) 

^ e-' In (^1 + Gu{dz)a{dU) < e-^^^h. 



II ^1/2 



□ 



Proposition 4.3. Let h e L{rK{a, (3)). Then the distribution of the K-valued random 
variable Y'^ is infinitely divisible and has Laplace transform 



,-/3((7)r- 



dMm(d[/) (4.24) 



exp(- 



In 1 



(e,c/> 



Gc/(dz)a(dC/). 



(4.25) 



where for a-a.e. f/, Gj/ is a Thorin measure measure on IR+ which is the image of 
Lebesgue measure on R+ under the change of variable s — > l3{U)/h{s, U). Moreover, the 
Levy measure ofY'^ is 



UY>^iE)= f r lE{rU)^^dra{dU), E e B{K) 
Jsf„ Jo r 



where 



ku{r) = / e-'-'Guidz). 



(4.26) 



(4.27) 



Proof Using the obvious analogue for the Laplace transform of the formulae for the 
characteristic functions of the integrals with respect to i.d.i.s.r.m.s in ||9l l30ll we obtain 



Lyh (6) — exp — 



K Jo 



exp 



Js'f,, Jo 



( 



Vf,^ {dx)ds 



1 -e" 



rh{s,U)(e,U) 



-rli(U) 



-dra{dU)ds 
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As in the last proposition, let Gu{dz) be the measure on which is the image of the 
Lebesgue measure on M+ under the change of variable s — > j3{U)/h{s, U) = z{U). Then 




Ly.(e) = exp ( y I I [l- e-'-^Q'^) j dra{dU)d 

(l - e-''<Q'^>) ^Gu{dz)draidU)^ 



Hence, combining 1 2.2 1 w^ith the existence conditions for the integral, Gjj[dt) is a Thorin 



measure on IR+ for a-a.e. U . □ 

4.4 Characterisation of Cone Valued GGC 

In this section we define Generalized Gamma Convolutions GGG(K) in the cone 
K and characterize this class as the distributions of the ii' -valued random elements 
represented by the stochastic integral ( |4.15[ l. The result is an extension to the cone val- 
ued case of the Wiener-Gamma integral representation of one-dimensional generalised 
Gamma convolutions, see Section 1.2 in 01611 . 

Similar to the multivariate case (see |[3]|), we define GGC(K) as follows 

Definition 4.4. The class GGC{K) is the collection of all infinitely divisible distribu- 
tions on K with Levy measure having a polar decomposition 

^m(^)=/ / lE{rU)^^draidU), E^B{K), (4.28) 

JS|IH JO 

where ku{r) is a measurable function in U and completely monotone in r for a-a.e. U. 



A probabilistic interpretation of the class GGC{K) is provided by Theorem 4.6 below. 
Proposition |4.3| says that the class of distributions of the Wiener-Gamma integrals 
Yh, h e L{rK{a,f3)), are GGCiK). 

We now prove that all distributions in GGCiK) have a Wiener-Gamma integral rep- 



resentation. For simplicity we consider the case without drift, that is 'I'o = in (4.1 
and (|4.3|l. Otherwise 



Theorem 4.5 (Wiener-Gamma characterization of GGC(K)}. For any fixed Borel-meas- 
urable function (3 : Sy^n M+ bounded away from zero it holds that 

Y'' ^ J h{s,x/\\x\\)xN{ds,dx) : a a finite measure on B{Sf,^^),h e L{rK{a, jS))^ 
= GGCo{K) 

with GGCq{K) denoting all generalized Gamma convolutions on K without drift. 

The condition on /? above is needed to ensure the existence of the Gamma random 
variables for all finite measures a. The result implies that starting with any fixed homo- 
geneous (or non-homogeneous with /3 bounded away from 0) Gamma distribution one 
can obtain all generalized Gamma convolutions as the sum of a constant and a Wiener- 
Ito integral with respect to the jump measure obtained from this fixed distribution. 
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Proof. Let e GGCq{K) with Levy measure given by (4.281. Since fc[/(r) is com 



pletely monotone in r for a-a.e. U, there exists a Radon measure Gu such that ku{r) 



Jo 



^Gu{dz). Moreover 



mm < oo. 

r 



(4.29) 



Let Fq^{x) — Guidz) for x > and i^Q^(s) be the right continuous generalised 
inverse of Fgu{s)- Let h{s, U) = l/i^c^(s) and h{s, U) = f3{U)h{s, U) for .s > 0. It follows 
as in the one dimensional case that Gu^, a-a.e. U, is a Thorin measure which is the 
image of Lebesgue measure on (0, oo) under the change of variable s — > l/h{s, U). That 
is. 



_£(£)x 

e ds = I e 

Jo 



'o 



'-Gu{dz), x>0 



and 



CO /'OO 



Jo 

oo poo 



r 

h{s,U) 



JO m) 

h{s,U) 



_ ^-Piu)r/Hs,u)^^ drdsaidU) 
dsa{dU). 



Thus ^^?29\ , Proposition [O] and (|4]20]l imply h e L(rif(a,^)). Let A^(ds,da;) be the 
Poisson random measure associated to F/f (a, (3) and Y'^ = h{s,x/ ||a;|l)a;A^(ds, dx). 

Then Proposition 4.3 shows that /x is the distribution of Y'^ which concludes the proof. 

□ 

We also have another characterization of GGG{K), similar to a characterization of 
multivariate GGC proved in ||3l Theorem F]. This gives another probabilistic interpreta- 
tion of GGG{K). 

We call XV an elementary Gamma variable in X if X is a non-random non-zero 
vector in K and y is a non-negative real random variable with one-dimensional Gamma 
distribution r(a, j3). 

Theorem 4.6. GGG(K) is the smallest class of distributions on K closed under convo- 
lution and weak convergence and containing the distributions of all elementary Gamma 
variables in K. 



Proof. The proof is along the same lines to that of Theorem F in fT, p. 27]. 



□ 



This implies also that GGC(X) is the smallest class of distributions closed under 
convolution and weak convergence containing all iiT -valued F distributions in the sense 
of Definition |4.1[ It is trivial to see that GGC(if ) includes all /v -valued stable distribu- 
tion using the spectral representation of stable Levy measures and maps h of the form 

h{s, U) = {se{U)T{a + l)}-^ with < a < 1. 



5 Positive Definite Matrix Gamma Distributions 

In this section we consider the important case of non-negative definite Gamma ran- 
dom matrices. This corresponds to the closed cone K = oi symmetric non-negative 
definite d x d matrices with inner product {X,Y) = tr(X^F), X,Y ^ S^. When X is 
in the open cone we write X > 0. When dealing with random matrices, a useful 
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1/2 

matrix norm is the trace norm defined for X e Md(M) as ||X|| = tT{{X^X) ' ). We 

write S+ I = sf^^ = S||.|| n S^. For X e S^J", \\X\\ = tr{X) and, in particular, if U e S+ |, 
tr(C/) = \\U\\ = 1. By Proposition 3.4 it is not important which norm we use. So we 
choose the one most convenient to work with. 



5.1 General Case 

The matrix Gamma distribution ^ ^ (a, f3) on §^ has the Laplace transform 

L^(e) = exp<^ - / / i^l - e-'-^'^^^^j dra(dC/) LveeSj (5.1) 

with alternative representation 

L^{e) = cxp I - / In (l + «(dt/) | ,Ve G §+. (5.2) 



+ 



Additional properties of Gamma random matrices to those for the general cone val- 
ued case in Section 4.2 are the following. 

If M is a symmetric random matrix with Gamma distribution rg+ (a, (3), tr(M) follows 
a one-dimensional Gamma convolution law. However, in the homogeneous case P{U) = 
(3o > 0, tr(7\f) has a one-dimensional Gamma distribution r(Q;(S||^jj), /3o). 

Proposition 5.1. a) If M ^ T^+{a,l3), tr(M) has a one-dimensional GGC law with 
Laplace transform 

Eg-etr(Af) =exp|-^°°ln(^l + ^^t;„,^(ds)| (5.3) 

where Va,i3 is the Thorin measure on (0, oo) induced by a{dU) on under the trans- 
formation U — > /?([/). 

b) If M ^ rg+ (a, /3) with j3{U) = /3o, tr(M) has the one-dimensional Gamma distribu- 
tion T{a{Sj A, Po). 



Proof For 6* > 0, let 9 = BU- Since Ee-^*'(^^) = L^(e), from (4.11 



= exp I - y In ^1 + Va^p (ds) 



where Va.j3 is the measure on (0, 00) induced by a(dJ7) on S||^|| under the transformation 
U — > P{U). Then, using (2.11 we obtain (a). For (b) we observe that from the first 

equality in the last expression with I3{U) = /3o, we obtain Ee"'"'^^^) = (1 + 6i//3o)""'^ii . 

□ 

Any matrix Gamma distribution rg+(a,/3) is self-decomposable and if supp(a) is of 
full dimension, it is absolutely continuous with respect to the Lebesgue measure on 
(which can be identified with M'*(''+i)/2) and so there is a density. The proof follows 
from the multivariate case, identifying the cone with Moreover, since the 

Lebesgue measure of S^\Sj is zero, the distribution T^+{a,P) is supported in the the 
open cone X > 0. In other words 
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Corollary 5.2. Let M he a random matrix with Gamma distrihution r^+fa, /3) with 
supp(a) of full dimension. Then V{M > 0) = 1. 

The following result is an adaptation of Proposition 3.8 to special linear operators 
preserving the cone Observe that all invertible surjective linear operators preserv- 
ing §+ are of the form X CXC'^ with some C e GLd{R) (see [20..21,l). 
Proposition 5.3. Let M ^ r=+(Q!,/3) with respect to the trace norm \\-\\ and let C e 

gCd{R). Then Y = CMC^ ~ rg+ (ac,/3c), where Sj,^^^ = S||.||^ n for = 
\\C~^BC-'^\\ and 

ac{E)^a{C-^EC-^), G S(S+ |^), (5.4) 

and 

l3c{U) = PiC-^EC-^), VC/eS+i^. (5.5) 
Moreover, Y ^ T^+(ac, f^c) with respect to the trace norm || || where 

Sc(i?)=/^^^^^li.(^)ac(dC/), yEeBiSl,) (5.6) 

and 

dc{U) ^ (3c ( J^) \\U\\c , VC/eS+|. (5.7) 



Italic. 

Example 5.4. (Diagonal matrix with independent entries). As pointed out in |[5|, an 

infinitely divisible non-negative definite random matrix M has independent components 
if and only if it is diagonal and therefore its Levy measure is concentrated in the diago- 
nal matrix axes i?*' e i = 1, ...d. Thus, a Gamma random matrix M ^ r^+{a,f3) has 
independent components, if and only if there exist non-negative numbers Pi, ...,13d such 
that the Levy measure vx is given by 



<i />oo —rf-ii 

um{E) = Va(-B..) / lE{rU)- dr E e S(§^). 



Further examples are considered in the next section. 

5.2 The Ar-distribution 

We now introduce a special matrix distribution ATd{rj, S) in the open cone §J with 
parameters, r; > (rf — 1) /2 and S S Sj. We study several properties including a relation 
between cumulants of AT{ri, E) and the moments of a Wishart distribution. 

The multivariate Gamma function, denoted by Td{ri), is defined for Re(?7) > {d— l)/2 

as 

Tdiv)= f e-*'-(^)|X|''-(''+^'/'dX, (5.8) 
Jx>o 

where dX is the Lebesgue measure on §J (identified with Rdid+i)/2y^ 

see for example 

0241 p. 61]. An alternative expression for Tdirj) is ( Ii24t Theorem 2.1.12]) 



Tdiv) = ^^("^-i^/^nr (v-li^- 1)) , Mv) >{d- l)/2. 



i=l 



(5.9) 



The special infinitely divisible matrix Gamma distribution ATdij], E) is defined as 
follows. For 77 > (d — l)/2, consider the measure p^(dX) = grj{X)(iX on the open cone 
Sj" where 

g,{X)=Cd^ \Xr^'+'^^\ X>0, (5.10) 

(tr(X)) 
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and 



(5.11) 



and LJd.rj > is given. 

Proposition 5.5. Let rj > {d— l)/2. There exists a homogeneous Gamma matrix distri- 
bution rg+(a^, /3) with respect to the trace norm where /3{U) = 1 for each U e S||;|j and 
a,, is the measure on Sf,, given by 



dU 



\U\ 



{d+l)/2 



(5.12) 



with a,,(S|f;j|) — ^d,i]- Moreover, the Levy measure of rg+(a,,,/?) is p,f and has a polar 
decomposition 



PviE) 



l-(-t^)'7^d.J^^, EeB{^:). (5.13) 



Proof. To show existence of the matrix distribution rg+(a^,/3), by Proposition 3.3 it 
suffices to prove that a,, is a finite measure, since trivially (3 satisfies (|4.12|l. The fact 
that rg+(a,,,/?) is concentrated in the open cone §J will follow by Corollary 5.2 since 



from 1 5.12 I supp(a^) has full dimension. 
For X > make the change of variable 



X = rU,r = tr{X), tr([/) = l,dX ^ r'^'-'^+^^/'^-^drdU 
e fact that \rU\ = 



(5.14) 



( Il23l p. 111]). Using this in \5.8\ and the fact that \rU\ ^r'^\U\ 

dU 

r"-''^e'' (lr\U[' — 
dU 

wv — 

Is 



\u\ 



{d+l)/2 



ic/r 



(5.15) 



and hence a^(S^||) = ujd,ri- Using again the change of variable (5.141 we have (5.13 1. □ 



Definition 5.6. Let ?/ > (d — l)/2 and E e Sj". An infinitely divisible p x p positive 
definite random matrix M is said to follow the distribution ATdir], S) if it has Gamma 
distribution rg+(a^,s,/^s) with respect to the trace norm where I3^{U) = tr(E~^J7) and 



lr,,s(dC/) = 



1 



|I]|"tr(E-iC/)'"' 



a„(dC/) 



(5.16) 



and is given by 1^5.12 ). 

Remaric 5.7. a) The distribution ATd{T],T,) has also as a parameter the total mass 
a,,(S^I|) = LJd.rj- This parameter is conjectured to be of particular importance when 
considering the limiting spectral (eigenvalue) distribution as the dimension goes to in- 
finity since it may then depend on r; or d. Particularly, interesting choices of w^.,, in this 
connection should be dr], d or a constant. 

b) The case r/ = (rf + l)/2 was considered in Barndorff-Nielsen and Perez-Abreu (Sj. 

c) It is interesting to note that for 77 g ((d— l)/2, {d + l)/2) the Levy density becomes 
infinity at the non-invertible elements of (i.e. the matrices which are positive semi- 
definite, but not strictly), whereas for 77 > {d+l)/2 the Levy density becomes zero at the 
non-invertible elements of For r/ = (d + l)/2 we have that a^^/^ is the uniform mea- 
sure on the unit sphere. This observation is also interesting in relation to applications 
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using e.g. ATd{rj, E) matrix subordinators. If all jumps should really be strictly positive 
definite , then a model with 77 > (rf+l)/2 should be appropriate, whereas one should use 
77 S ((d— 1)/2, (d+l)/2), if it seems desirable to have most of the strictly positive definite 
jumps rather close to non-invertible matrices. The latter should be especially useful in 
stochastic volatility models in finance where very often news (resulting in jumps of the 
covariance) should only really affect a single asset or a special group of assets (like 
stocks of companies from the same country or the same branch of industry). 

Note that if M - ATi-q, U), then E^/^Af i;i/2 _ ATair], E). This follows from Proposi- 



tion 5.3 which also gives together with Proposition 5.5 that ATdir], E) has Levy measure 



p^,s(dX) = g,^^s{X)dX where 



The existence of moments of all orders of AT\i{r],lj_) follows since l |4.12 1 is trivially 
satisfied. The same is true for ATdir], E) since T.^/'^MT.^/^ - ATdiv, E). 

In the homogeneous case the distribution AT^ir], aid) with a e M+ is invariant under 
orthogonal conjugations and the trace follows a one-dimensional Gamma distribution. 

Lemma 5.8. Let 77 > (rf - l)/2 and ct > 0. 

The distribution ATdir], aid) is invariant under orthogonal conjugations, 
b) If M ~ ATd{ri,ald), then tr(A/) follows a one-dimensional Gamma distribution 

T{LOd,r,,Cr). 

Proof. It is well known that the measure dX/ \x\'^•'^^^^^^ is invariant under the conju- 
gation X — > CXC^ , for X > and any non-singular matrix C (see [jL3, Example 
6.19]). The determinant and the trace norm functions are invariant under the conju- 
gation X 0x0^ , for X > and any O e 0{d). Thus the Levy measure Pr^^aU with 
density ( |5.17| l, E = aid, is invariant under orthogonal conjugations and so the matrix 
distribution Ar{r],ald) is. Proposition |5 . l[ h) gives (b). □ 

The cumulants of the distribution Ard{rj,ald) are related to the moments of the 
Wishart distribution, as we prove below. Recall that a d x d positive definite random 
matrix W is said to have a Wishart distribution Wd{n, d) with parameters n > d — I and 
E e Sj', if its density function is given by 

fw{A) = ^e-^t^s-^A) |^|(n-<i-i)/2 ^ Q (5_^g) 

2*r,(f)|E|^ 

As usual we denote hy Ai^B the tensor product of the matrices A and B. We recall that 

if A and B are in M„, then tr(yl (g) B) ^ tr(A)tr(B) and \ A(g)B\^ | \B\" . 

We use the notation B{n,m.) — r(n)r(r7i)/r(n -I- m) for the Beta function, where 
Re(n) > 0, Re(TO) > 0, and for the multivariate Beta function, denoted by Bd{x, y), 

Bd{n,m) ^ ^din)rd{m) j^g(^) > q, Rc(?7i) > 0. (5.19) 
Td(n + m) 

Proposition 5.9. Let rj > (d - l)/2 and i>e the Levy density of ATdii], S). Let 

W be a random matrix with Wishart distribution Wd{2ri, E). For any integer p > the 
following three identities hold: 
a) 

f XPg^MX)dX = ^B{ijd,p)EWP. (5.20) 
Jx>o ^ 
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b) 



c) 



X>() 



^,,s(X)dX = ^BM,p)E(H^«f). 



X>0 



\Xfgr,MX)dX 



2pd 



^B{rjd,pd)E{\Wn 



B{r]d,pd) p 



(5.21) 

(5.22) 
(5.23) 



Proof. The existence of the integrals l |5.2Q[ l-l [5.22[ l is seen as follows. The finiteness of 
the p-th moment of ATd{ri, S) is equivalent to the existence of the p-th moment of its 
Levy measure (away from the origin). Then /^^q \\X\f grj,^{X)dX < oo for any p > 
gives the existence of the integral in \5.20\ . Since for X > tr(X®P) = ~ (tT(X )f ^VQf 
and \X\ < \\X\f, one also obtains the existence of the integrals in l5.21l and 1 5.22 1, 
respectively 

The identities in (a)-(c) are consequences of a more general result for g-homo- 
geneous functions which we now prove: Let h : ^ H he a function such that 
h{rX) = r'^h{X) for any r > and X e and some fixed g > 0, where H is §J, 
(Sj)*^^ and (0, oo) for (a),(b) and (c) respectively. In the cases (a) and (b) we have q = p 
while for (c) q = dp. 

The change of variable V = E^^/^XE^^/-^, the invariance of the measure 



\X\ 



<d+l)/2 



dX under non-singular linear transformations and writing h^{Y) 



/l(El/2rEl/2) give: 

hiX)g^^^iX)dX = Cd,r, 



J = 



X>0 



h{X) 



-tr(S-^X) 



g-tr(y) 



x>o 1^ 

V-{d+l)/2^y^ 



(tr(S-iX)) 



rjd 



\x 



,n-{d+i)/2 



dX 



Using l |5.14[ l, the definition of the Gamma function, first for T{q) and then for T{T]d + q), 
and ( |5.1ip give 



Cd,r,r(g) 

cd^v^jq) 

uJd,v^ivd,q) 

^div) 
uJd,r,^ivd,q) 

2v'iTd{r,) 
uJd,n^ir]d,q) 



/is(c/)it/r 



dU 



|JJ|(d+l)/2 

r'"^+'?-i/is(C/)e-'^dr |f/r 



dU 



|f/|(rf+l)/2 



V>0 



-A 



-itr(A) 



\A\' 



dA 



(d+l)/2 



A] fw{A)dA = LOd^r^'Qivd^q) 



(5.24) 
(5.25) 
(5.26) 



where in l |5.24[ l we used again l 5.14 l, in (5.251 the change of variable ^4 = V/2 (with 
dV = (i/2)('^+i)/2dA), and for ( |5.26| l the fact that fw is the density ( |5.18| of the random 
matrix W with Wishart distribution Wd{2ri,Y,). 

Then (a), (b) and the first equality in l |5.22[ l are proved. The second equality in 
( |5.22[ ) follows using the fact that E \W'p\ = d{v + p) d{v) (see Muirhead Ell p. 

101]). □ 
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In particular, the mean E{M) and covariance Cov(M) = E(M® M) -E(M)®E(Af) of 
M ^ ATdirj, S) are expressed in terms of the mean E(Vl^) and the second tensor moment 
E{W)<SiE{W) of the Wishart distribution Wd{2r], E). Recall that the commutation x (f 
matrix K is defined as 

d 

where Hij denotes the c? x d matrix w^ith hij = 1 and all other elements zero. The m-th 
moments and cumulants of a d x d random matrix are d^^-dimensional objects which 
need to be represented in a concise and at the same time easy to handle way. As usual 
for random matrices we define the moments and cumulants using the tensor product, 
e.g. the m-th moment of a random matrix X is understood to be E{X®"^). An alternative 
would be to use the vec-operator to transfer the matrix into an element of first, but 
typically this leads to formulae that are more cumbersome to handle. Now we have: 

Corollary 5.10. The cumulants of the random matrix M ~ ATdirj, S) are proportional 
to the tensor moments of the Wishart distribution. In particular 

E(Af) (5.27) 

and the matrix of covariances between elements of M is given by 

Cov(M) . ....^^^ ((l + ^) Id^ + i^) (S « S). (5.28) 

Proof. The first assertion follows from (b) in Proposition |5.9[ Since the first moment of 
M equals its first cumulant and its matrix of covariances equals its second cumulant. 



then (5.271 follows from (a) in Proposition 5.9 with p = I and since E(VF) — 27]!^ for 



W - VFd(2r;, S). From [24, p. 90] we have 

Cov{W) = 277(I<j2 +K){Y,(^Y.). 
Using (b) in Proposition |5 . 9| with p = 2 we have 



Cov(M) = [ X®2g^,.s(X)dX 
Jx>o 



^ = E(W^^) 

4 {nd+l)r]d ^ ' 



^d,r] 1 

4 (nd + l)r]d 



Cov{W) + [E{W)f^ . 

Hence l |5.28| l follows. □ 



In particular, when w^.,, = drj, E{M) = rjT., as in the Wishart case. On the other hand, 
when ujd,rj = d, E{M) = E. 

This result is of particular importance in applications, since it implies that the second 
order moment structure is explicitly known which may allow method of moments based 
estimation of models using ATdir], E) matrix subordinators as the stochastic input (e.g. 
II291I1 

The following result states an interesting relation with the so-called Marchenko- 
Pastur distribution of parameter A > 0. Recall that the moments of this distribution are 
given by (see 
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Lemma 5.11. Let e e R. For any integer p > Q, as d ^ oo and d/rj A > 



where Kp{\) = r{p)^p{2\). In particular, for e ~ 1 



lim 



Proof. IfW Wd{2ri, U), then 









-/ 


tr 




d Jx>o 







9v.iAX)dX 



\, ifp = 1 
0, ifp > 2 ■ 



lim -Etr 

d^oo d 



-Mp(2A). 



(5.30) 



(5.31) 



(5.32) 



as 277/d — J- 2 A due to the well known Marchenko-Pastur Theorem [2], for any p > 0, 
By the Stirling approximation T{z + 1) ~ V2TTz{z/e)^ for z ^ 00, for r/ and d large 



T{r,d + p) 



{nd)~P. 



(5.33) 



Using rj/d^ \ and ( 5.32 1 in ( 5.20 1 gives: 



x>o V d" 



9,jAX)dX 



0.,,, rM)r(p) 1 wp 
d'^p ri-qd + p) d ^ 2 ' 



T{p)ujdAvd)-''d-P' 



2 



-Etr 
d 



T{p)LOd,„d-P-P' 
Kp{\)LOd,nd^P^^+^\ for ?7, d large. 



2 



which proves the lemma. 



(5.34) 
(5.35) 

□ 



Conjecture 5.12. We conjecture that the above Lemma is a first step to study the 
asymptotic spectral distribution of the random matrix M ^ ATd{ri, S). More specifically, 
the right hand side of (5.3ip must Jbe related to the pth-cumulant of the pth-moment of 
the mean spectral distribution of M, which in turn should allow the identification of the 
limiting spectral distribution. 

5.3 Further examples 
5.3.1 BF-distributions 

Let d> I and q — 1, d be fixed. Consider the Levy measure on §^ given by 



where (3o > and 



,{dX) = 



ad,q{E) 



-0o\\X\\ 



X 



-ad,,{dX), XeSj\{Q} 



lE{rU)drad,g{dU). 



(5.36) 



(5.37) 



Here ad.q/dis the probability measure on the sphere Sj|;|| induced by the transformation 
V ^ U ~ , where the d x q matrix V is uniformly distributed on the unit sphere 
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of the linear space M^x^ClR) of dxq matrices with real entries, with the Frobenius norm 

WYWl^triY^Y). 

An infinitely divisible d x d symmetric random matrix M with Levy measure I'q has 
the Gamma distribution T^+{ad,q, (3q), since Vq has a polar decomposition 

p poo — r 

i^qiE)^ / lE{rU)—drad,q{dU), E e B{Ea). 

Jst „ Jo T 



We call this distribution the BTd{q, Po) distribution. 

Remark 5.13. a) We observe that the support of Vq is concentrated in matrices of rank 
q in S([. Hence this support is of full dimension. Then, by Corollary 5.2 T^+{ad,q, Po) 

has support in the open cone S^^. 

b) The case q = I was considered in Perez-Abreu and Sakuma II27II in the context of 
random matrix models for free generalised Gamma convolutions. They considered the 
Hermitian case for which working in the setup of Mdxq(C) is needed, but otherwise 
the above steps can be carried out in a straightforward way. 

c) For E e §t one can consider invertible linear transformations of T^+{ad,q, Po) 

d 

to obtain infinitely divisible positive definite matrix Gamma distributions rg+ (a^^s, /9s) 
with Levy measures of the form 

-||s-ix|| 

similar to the family F of matrix Gamma distributions considered in the last section. 

The following properties are easily proved. 
Proposition 5.14. Let M ^ r^+{ad,q, Po) and g = 1, ...,d be fixed. Then, 

a) M has an invariant distribution under orthogonal conjugations. 

b) E \\Mf < oo for any k > 0. 

d) tr(Af) has a one-dimensional Gamma distribution r{d,(3o). 
5.3.2 Matrix Gamma-Normal distribution 

In the one-dimensional case, the so called variance gamma distribution is popular in 
applications in finance, see II22I1 . This distribution is a mixture of Gaussians having 
a random variance following the one-dimensional Gamma distribution. As an appli- 
cation of the matrix Gamma distribution, we now present a matrix extension of the 
one-dimensional variance Gamma distribution. 

Let Z he a d X q random matrix with independent standard Gaussian distributed 
entries, i.e. 

Eexp(itr(e^Z)) = exp(-itr(e^e)), VG G M^xgW- 

Let X be a random matrix with the Gamma distribution rg+ {a, (3) and independent of Z. 

Consider the random linear transformation Y ~ X^I'^Z. Using a standard conditional 
argument we compute the characteristic function of the dxq matrix as follows: 

Ecxp(rtr(e^y)) = ExEz cxp(itr(e^X^/2z)) X) 

= E;,|exp(-itr(xV2ee^Xi/2))| 



= E^ |cxp(-^tr(ee^X))| 
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Then, using ( |4.11[ l 

Eexp(ztr(eTy)) = exp|-^^ In (^1 + l^^i^^^^ a(dt/)| , (5.38) 

for each e Mdxg(]R). Using the terminology in |[6]|, we can say that Y has a MatG 
distribution, which is infinitely divisible in Mdxq(]R). 

Similar to the one-dimensional case, we call this distribution the matrix Gamma- 
Normal distribution with parameters a and (3 or more specifically the d x g-dimensional 
matrix rg+ {a, /3)-Normal distribution. We observe that Y has a symmetric distribution 

in the sense that — Y^'^Y and also that Y has a distribution invariant under orthogonal 
conjugations if I3{U) = /?o and a{dU) is invariant under orthogonal conjugations. 

Remark 5.15. If P{U) = /3o > and q — d, tr(Y') has a one-dimensional variance Gamma 
distribution with the following characteristic function: for 6* e M, O = Old, 



EexpriMrfF)) = exp i - / 



In ( 1 + ) «(dC/) 



1 



1 

1 



a2 



-q(S 



+ 



"(S^.„) / 1 \ -"(St,,) 



2/3o J V ^/Wn 

Thus, tr(F) has the same distribution as — V, where V has a one-dimensional Gamma 
distribution r(a(S||^||), ^/2/3o) and V^'is an independent copy of V. 

Remark 5.16. Let q ^ I, i.e. the resulting matrix Gamma-normal distribution is M''- 
valued. Thus for 9 e W'-, QQ^ has rank one. 

Assume additionally that the measure a is concentrated on the rank one matrices, 
that is U = uu^ with u e R"^ (and the first non-zero component of u being positive, to 
make u unique given U). Let a be the measure on the unit sphere Sgd ||.|| of M'* induced 
by a under this transformation. Using this we write the integral in the right hand side 
of ( |5.38| l as follows 

, / itr(c/ee^)\ , 




(5.39) 

Interestingly, l |5.39| l implies that the matrix Gamma-normal random variable can be 
represented (in this special case) as Xi — X2 with Xi,X2 ^ Td{a, /3) being independent 
where /3 = a/2^(uuT). 

Hence, the matrix Gamma-normal distribution with q ~ \, which can indeed be 
regarded as a d-dimensional generalisation of the univariate variance Gamma distribu- 
tion, inherits interesting properties well-known in the univariate case. 
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